Beyond Carbon K-edge harmonic emission using spatial and temporal synthesized 

laser field 
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We present numerical simulations of high-order harmonic generation in helium using a tempo- 
rally synthesised and spatially non-homogeneous strong laser field. The combination of temporal 
and spatial laser field synthesis results in a dramatic cut-off extension far beyond the usual semi- 
classical limit. Our predictions are based on the convergence of three complementary approaches: 
resolution of the three dimensional Time Dependent Schrodinger Equation, time-frequency analysis 
of the resulting dipole moment and classical trajectories extraction. Employing a combination of 
temporally and spatially synthesised laser field provides coherent XUV photons beyond the car- 
bon K-edge which is of high interest for initiating inner-shell dynamics and study time- resolved 
intra-molecular attosecond spectroscopy. 
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Near-edge X-ray absorption spectroscopy [H, E| is a 
very powerful technique to probe local chemical envi- 
ronment of molecules. Providing to such technique an 
attosecond time resolution will open a new route to ex- 
plore inner shell inducing ultra-fast dynamics of charges 
in molecular system and therefore it is expected to be 
an important breakthrough. One way to achieve attosec- 
ond X-ray is the generation of high harmonics (HHG) 
driven by mid-infrared strong laser field, as it is well 
known that the HHG cut-off scales as A 2 (A being the 
laser wavelength) [3|. A first demonstration has been 
very recently done |4| where high requirement on the 
mid-infrared wavelength femtosecond laser source may 
remain quite technologically challenging. In addition, 
it has been demonstrated that the harmonics efficiency 
strongly decreases with longer wavelength according to 
a law in A -5 ' 5 0-0] ■ Therefore high demands are on 
developing alternative routes based on existing conven- 
tional femtosecond laser sources at 800 nm. So far a large 
amount of strategies to tailor such input laser pulses for 
HHG have been carried out. Using a combination of 
two or more colours laser fields to synthesise temporally 
the final strong laser field jl-18|; or using chirping tech- 
niques . The new approach we propose consists in the 
combination of two techniques to shape the final laser 
field both in time and in space. The temporal synthe- 
sis is obtained using two few-cycle laser pulses delayed 
in time [20| . Few-cycle laser pulses are now commonly 
used and are produced using post-compression stage for 



22[. The spatial synthe- 



attoscience applications |21l . 
sis is obtained by using a non-homogeneous laser field 
as produced by nano-plasmonic antennas for in- 
stance. The goal of this paper is to show that synthesising 
the laser field both in time and in space is a new route 
to generate photon via HHG beyond the carbon K ab- 
sorption edge and this using near infrared wavelength as 



800 nm, which corresponds to conventional femtosecond 
Ti:Sa laser sources for strong field physics. 

To synthesised temporally the laser field we use a 
method described in [20j where two delayed pulses of 4- 
cycles at 800 nm are combined and correspond to, 
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where E is the laser electric field amplitude, in atomic 
units (E = y/T/To with I = 3.5 x 10 16 W/cm 2 ), w = 
0.057 a.u. the laser frequency corresponding to A = 800 
nm, T the laser period, N the total number of cycles in 
the pulse, <fi the carrier-envelope phase (CEP) and r the 
time delay between the pulses. In our simulations we 
consider the case N = 4 and = 0. 

According to our numerical simulations the optimal 
time delay between these two pulses corresponds to r = 
1.29T. This way the laser amplitude after the superpo- 
sition is equal to the two input pulse replicas. The two 
replica have the same carrier envelop phase (CEP), which 
is always the case if the replicas are produced from the 
same initially CEP stabilised source. With such config- 
uration an extension of the cut-off has been predicted 
up to 4.5£/p, where U p is the quiver energy associated 
to the laser field employed (U p — xj^)- Most of the 
numerical and semi-classical approaches to study HHG 
in atoms and molecules consider ho mog eneous temporal 
and spatial laser field distributions [26|, |27J . Some stud- 
ies have study the effect of homogeneous temporal distri- 
bution together with synthesised spatial profile. In [28| 
they proposed to synthesised a spatially flat top profile 
to produce isolated attosecond pulses showing that the 
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FIG. 1. 3D TDSE harmonic spectrum in Helium atom gen- 
erated by laser pulses described in Eq. @ and P = 0.002 for 
a plasmonic enhanced intensity I — 1.4 x 10 15 W/cm 2 . 



spatial domain control of the laser field provides degree 
of freedom to control HHG. 

Our idea is therefore to use not only temporally synthe- 
sised laser field using a combination of two CEP stabilised 
few-cycle pulses as described by Eq. ((T|) and © (referred 
to as double pulse configuration) but also spatially lo- 
cally synthesised (referred to as non-homogeneous con- 
figuration). This non- homogeneous spatial distribution 
of the laser field can be obtained exp erimentally by us- 
ing a non- homogeneous laser field 23-25J as produce d by 
nano-plasmonic antennas |29l| . m etallic waveguides [30|, 



metal [31. l32j a nd dielectric |33j nanoparticles or metal 
nanotips The laser electric field is then no longer 

homogeneous in the region where the electron dynamics 
takes place. Very recently, studies about how HHG spec- 
tra are modified due to non-homogeneous fields, such as 
those present in the vicinity of a nanostructure irradiated 
by a short laser pulse, have been published 23-25, 
leading to the so called 'plasmonic enhanced intensity' 
showing the increase interests such new field rises . 

In order to calculate the harmonic spectra we re- 
solve the three dimensional Time Dependent Schrodingcr 
Equation (3D-TDSE) in the length gauge while employed 
a double pulse and non-homogeneous driving laser field. 
The harmonic yield of an atom is proportional to the 
Fourier transform of the dipole acceleration of its active 
electron and it can be calculated from the time propa- 
gated electronic wave function. We have used our code 
which is based on an expansion in spherical harmonics, 
YJ m , considering only the m = terms due to the cylin- 
drical symmetry of the problem. The numerical tech- 
nique is based on a Crank-Nicolson method implemented 
on a splitting of the time-evolution operator that pre- 
serves the norm of the wave function. We base our stud- 
ies in He due to the fact that a majority of experiments 
in HHG are carried out in noble gases. Hence we have 
considered in our 3D-TDSE code the atomic potential 
reported in [i(J to describe accurately the He atom. We 
remark that our scheme can be applied to any atom or 
molecule, once the adequate potential is chosen. The 
coupling between the atom and the laser pulse, linearly 
polarized along the z axis, is modified in order to treat 



spatially non-homogeneous fields as follow: 

Vi(z,t,r) =E{z,t,T)z (3) 

with 

E(z,t,r) = E(t,T)(l+pz) (4) 

where E(t, t) is the laser field defined by the sum of the 
pulses defined in Eq. (JTJ) and (|2} , the parameter j3 defines 
the strength of the non-homogeneity and the dipole ap- 
proximation is preserved because /? -C 1. The harmonic 
spectrum then obtained in He atom for a j3 — 0.002 is 
presented in Fig. [T] We can observe a considerable cut- 
off extension up to 12.5U P which is much greater than 
whilst the double pulse configuration is employed alone 
(that leads to a maximum of 4.5C/ p [20[ ). This large 
extension of the cut-off is therefore a signature of the 
combined effect of the temporal double pulse and the 
spatial non-homogeneous character of the laser electric 
field. For this particular value of laser peak intensity 
(1.4 x 10 15 W/cm 2 ) the highest photon energy is greater 
than 1 keV. Note that the cited intensity is actually the 
plasmonic enhanced intensity, not the input laser inten- 
sity. The latter could be several orders of magnitude 
smaller, according to the plasmonic enhancement factor 
(see e.g. j^HHU). 

In order to understand the mechanism link to this par- 
ticular extension of the cut-off, we investigated the clas- 
sically retrieved trajectories by solving the Newton equa- 
tion of an electron in the such laser field Eq. We cal- 
culated the energy of the returning electron as a function 
of the time duration of the laser pulse 25|. In Fig. [2] we 
show the classical rescattering energies (in eV) as a func- 
tion of ionization time (in blue) and recombination times 
(in black) including the effect of the non-homogeneous 
field varying the j3 parameter. The direct effect is then 
that the number of recombinations decreases with (3 in- 
creases. The non-homogeneity of the laser field acts 
therefore as an electronics trajectories filter, with the ad- 
vantage that the most energetic trajectories are selected 
leading to a cut-off extension greater than 1 keV. Fur- 
ther to this selection and extension of photon energies, 
an additional very interesting effect linked to the non- 
homogeneity of the laser field rises. Indeed for the case 
of ft = 0.002 the so-called "short" and "long" classical 
trajectories 41-43J recombines now almost simultane- 
ously. In other words, the non-homogeneity of the laser 
field acts as a temporal lens which forces to the electrons 
previously ionized at different times to recombine at the 
same time as it is shown in the panel (d) of Fig. [2] This 
way the attosecond chirp is significantly reduced. 

In addition, in Fig. [3] , the recollision time t r of the 
electron is presented as a function of the ionization time 
ti and for several values of /3. By considering ionization 
times of 1.25 < t, < 2.25 the long trajectories are those 
with recollision times t r > 2.5 optical cycles and, only 
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FIG. 2. Rescattering energies of electrons as a function of 
the ionization time (blue) and recombination time (black), 
for laser pulses described in Eq. Q for different values of 
/3 = 0.0 (a), p = 0.001(b), P = 0.0015(c) and p = 0.002(d). 
In this case the plasmonic enhanced intensity is / = 1.4 x 10 15 
W/cm 2 which corresponds to the saturation intensity for He. 
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FIG. 3. Recombination times as a function of ionization 
times for three of the four cases exposed in Fig. [2] for the 
He atom. The case p = is plotted in grey, P = 0.001 in blue 
and P = 0.002 in red. 



for the homogeneous laser field case (grey lines in fig- 
ure), these trajectories are clearly visible. On the other 
hand, short trajectories are characterised by t r > 2.5 
optical cycles and these are present for both the homoge- 
neous and non-homogeneous cases. The resu lts p resented 
in here are consistent with those shown in [24L 251, but 



note that now the laser pulse has a more complex tem- 
poral shape. These results on the fact that the long tra- 
jectories are modified both by the non-homogeneity and 
the temporal double pulse configuration, namely, the ho- 
mogeneous long trajectories (grey lines) with ionization 
times ti around 1.25 and 1.75 optical cycles merge into 
unique trajectories. The branch with ti ~ 1.75 now has 
ionization times more than half an optical cycle that get 
smaller whilst /3 increases. As a result, the time spent 
by the electron in the continuum increases. The electric 
field strength at the ionization time for short trajectories 
being greater than for long trajectories, and considering 
that the ionization rate is a non-linear function of this 
electric field, long trajectories are then less efficient than 
the short ones. On the other hand, short trajectories 
are almost independent of /3 and get noticeable different 
only for really high values of (3. Similar behaviour can 
be found in Fig. [3] for ionization times of 2.5 < ti < 3.5. 

This effect has important consequences in the har- 
monic spectrum plotted in Fig. [T] which exhibits a wide 
zone where no clear harmonics peaks are visible resulting 
in a continuum. As demonstrated in |21[ , for continuum 
obtained around 90 eV using a 800 nm laser field using 
single few-cycle pulse and homogeneous field, continuum 
generation in the harmonics spectra is a necessary con- 
dition for the generation of an isolated attosecond pulse. 
So far no clear method has been proposed to allow the 
generation of isolated attosecond pulse in the keV region. 
Therefore generating harmonics using the combine effect 
of temporal and spatial synthesis is a completely new 
route towards the first generation of isolated attosecond 
pulses beyond the Carbon K-edge. 



As a final test to confirm the underlying physics high- 
lighted by the classical trajectories analysis, we retrieved 
the time-frequency distribution of the calculated dipole 
(from the 3D-TDSE) corresponding to the case of a non- 
homogeneous laser field using a wavelet analysis. The 
result is presented in Fig. [4] where we have superim- 
posed the classical rescattering energies (in brown) to 
show how the two calculations match. Analogously, in 
Fig. @] we show the time- frequency analysis in the case of 
p = 0.002 that corresponds to the spectra presented in 
Fig. [TJ The consistency of the classical calculations with 
the full quantum approach is clear and confirms that the 
mechanism of the generation of this 12.5U P cut-off ex- 
tension exhibiting a nice continuum is the consequence of 
trajectory selection and temporal lens effect on the recol- 
lision time whilst employing the combination of temporal 
and spatial synthesised laser field. 
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FIG. 4. Time-frequency analysis obtained from the full inte- 
gration of the 3D TDSE and superimposed (in brown) classi- 
cal rescattering energies for the case of P = 0.002 and a plas- 
monic enhanced intensity / = 1.4 x 10 15 W/cm 2 correspond- 
ing to the case of Fig. [5] 

In conclusion, using a laser fields with a temporal su- 
perposition of two identical few-cycle pulses delayed in 
time together with a spatial non-homogeneity, we demon- 
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strate that temporally and spatially synthcsiscd laser 
field can lead to a new route for night harmonic genera- 
tion, attosecond pulse production and in a further extend 
to strong field physics. We demonstrate that the main ef- 
fect of this synthesised laser field on the harmonics results 
in a considerable extension of the cut-off energy up to 
12.5£/p which present a nice continuum generation com- 
patible with the generation of isolated attosecond pulse 
beyond the Carbon K-egde. This effect is understood by 
analysing the trajectories involved in the process using a 
classical and a full quantum approaches. Both analysis 
converge to the same conclusion: trajectories are highly 
selected whilst using a laser field that consists of a com- 
bination of the double pulse temporal synthesis and the 
spatially non- homogeneity. And a 'temporal lens effect' 
shrink the recombination time window, resulting in a nice 
continuum generation in the extended harmonics spectral 
region. This new approach provides a unique route to the 
generation of coherent attosecond light source beyond the 
Carbon K-edge directly from input 800 nm laser system 
and therefore is of interest of strong field applications 
time-resolved X-ray spectroscopy of molecular system. 
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